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ABSTRACT 



We present the spectroscopy of red giant stars in the dwarf spheroidal galaxy Leo I, aimed at further constraining its chemical enrich- 
ment history. Intermediate-resolution spectroscopy in the Can triplet spectral region was obtained for 54 stars in Leo I using FORS2 
at the ESO Very Large Telescope. The equivalent widths of Can triplet lines were used to derive the metallicities of the target stars on 
the [Fe/H] scale of Carretta & Gratton, as well as on a scale tied to the global metal abundance, [M/H]. The metallicity distribution 
function for red giant branch (RGB) stars in Leo I is confirmed to be very narrow, with mean value [M/H]^ -1.2 and dispersion 
f[M/H] — 0.2. By evaluating all contributions to the measurement error, we provide a constraint to the intrinsic metallicity dispersion, 
°"[m/hj,o = 0.08. We find a few metal-poor stars (whose metallicity values depend on the adopted extrapolation of the existing calibra- 
tions), but in no case are stars more metal-poor than [Fe/H] = -2.6. Our measurements provide a hint of a shallow metallicity gradient 
of -0.27 dex Kpc~' among Leo I red giants. The gradient disappears if our data are combined with previous spectroscopic datasets 
in the literature, so that any firm conclusions about its presence must await new data, particularly in the outer regions. By combining 
the metallicities of the target stars with their photometric data, we provide age estimates and an age-metallicity relation for a subset 
of red giant stars in Leo I. Our age estimates indicate a rapid initial enrichment, a slowly rising metal abundance -consistent with the 
narrowness of the metallicity distribution- and an increase of ~ 0.2 dex in the last few Gyr. The estimated ages also suggest a radial 
age gradient in the RGB stellar populations, which agrees with the conclusions of a parallel study of asymptotic giant branch stars in 
Leo I from near-infrared photometry. Together, these studies provide the first evidence of stellar population gradients in Leo I. 

Key words. Galaxies: dwarf spheroidal - Galaxies: individual (Leo I) - Stars: abundances - Local Group - Galaxies: stellar content 



1. Introduction 

A first-order estimate of the distribution of stellar metallicities in 
resolved galaxies can be obtained from photometry of red giant 
branch (RGB) stars. However, young metal-rich stars have the 
same colours as older metal-poor stars, a phenomenon known 
as "age-metallicity degeneracy". Although improved photomet- 
ric metallicity estimates can be obtained by combining optical 
and near-infrared photometry, with reduced degeneracy effects 
dGullieuszik et al.l I2007L I2008I) . the problem is never entirely 
overcome, and age and metallicity inextricably contribute to the 
observed colour. 

More direct measurements of the stellar metallicity distri- 
bution function (MDF) are obtained from spectroscopy. The 
most accurate determinations come from high-resolution abun- 
dance analysis which provides information on the relative abun- 
dances of chemical elements, and therefor e the strongest con- 
straints to the galaxy evolution models (e.g. lTolstov et al.ll2003l 
lLanfranchi & Matteuccil I2007I) . High-resolution spectroscopy, 
however, is limited to the brightest stars in nearby dwarf 
galaxies, and becomes unfeasible for more distant galaxies 
even within the Local Group and using lOm-class telescopes. 



A viable alternative for obtaining metallicities for a large 
number of stars (as needed to derive statistically significant 
stellar MDFs) is using low- or intermediate-resolution spec- 
troscopy. The infrared Ca triplet (CaT) method, originally de- 
vised Jojneasvrrejnejallicitie^ 

ters ( Armandroff & Zinnlll988t rArmandroff & Da Costal 1 19911 
Rutledge et al.ll997bl) . has now been applied quite extensively to 
RGB stars in dwarf galaxies. Using this method, va lues of [Fe/H] 
for dwarf galaxies have been obtained for F ornax dTolstov et al 



200l | : lPont e t al.l2004lBattagria et al.l2 006l Carina dKoch et al 



2006y^lptor jjTolstov et al.ll2001l), LeoII dBosleret : alj|2007t 



iKoch et all l2007al). NG C 6822 dTolstov et al.l 1200 ll) and the 



LMC dCole et al.ll2005l) . 



Send offprint requests to: M. Gullieuszik 

* Based on data collected at the European Southern Observatory, 
Paranal, Chile, Proposals No. 69.D-0455 and 71.D-0219 



In this paper we present new CaT spectroscopy for the dwarf 
spheroidal (dSph) galaxy Leo I. Leo I is known to have formed 
the bu lk of its stars at an intermediate epoch (e.g.. lGallart et al.l 
119991) . Along with LeoII, Leo I is one the most distant dSph 
satellites of the Milky Way, for which the influence of tidal in- 
teraction with our Galaxy on evolution must have been more 
limited than for nearby galaxies such as Carina or Fornax. As 
such, deriving its chemical enrichment history is of foremost im- 
portance for our kn owledge of the evolution of dSph galaxies. 
iGallart et al.l d 19991) derived the metallicity and reconstructed 
the star-formation history (SFH) of Leo I from HST/WFPC2 
data - they found a metallicity ranging from [Fe/H] = - 1 .4 to 
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Fig.l. Digitized Sky Survey image of the Leol field, centred at 10''08"'28 s :i, +12°18'23" (J2000). Starred symbols mark the 4 
metal-poor member stars, the circles are other member stars, squares represent interlopers. 



[Fe/H] = -2.3. Using the same data set, iDolphin] d2002l) found 
higher metallicity, ranging fro m [Fe/H] = —0-8 to [Fe/H] = -1.2. 
From the colour of RGB stars. iHeld et al. d2000t) deri ved a mean 
metallicity [Fe/H] 1.6 on the lZinn & West! d 19841) scale. 

The metallicity distribu tion of RGB stars in L eol has re- 
cently been investigated by Bosler etaT] d2007l) and iKoch et all 



d2007bl) . iBosler et al] d2007l) used the CaT method to analyse 
Keck-LRIS spectra of 102 RGB stars, and found [Fe /H] = 
-1.34 on the [Fe/H] scale of ICarretta & Grattonl d 1997b (here- 
after, CG97). The authors also proposed a new calibration based 
on Ca abundance, yielding a mean metallicity [Ca/H] = -1.34 
(cr = 0.21) . Also u sing measurements of the Can triplet lines, 
iKochet alid2007bl) found a mean metallicity [Fe/H] = -1.31 on 
the CG97 scale for 58 red giants. In both studies, the MDF is 
well described by a Gaussian function with a lcr width of 0.25 
dex, and a full range in [Fe/H] of approximately 1 dex. 

Leol is an interesting target also because, to date, it is one 
of the few Local Group dwa rf galaxies showi ng scarce evi- 
dence of a population gradient. IHeld et al] d2000) found that the 
old horizontal branch stars of Leo I are ra dially distributed as 
the intermediate-age helium-burning stars; K och et al] d2007bl) 
found no significant metallicity radial gradient. Different conclu- 
sions were fo und in our compa nion paper, based on near-infrared 
photometry dHeld et al] 120091) . showing that intermediate-age 
asymptotic giant branch (AGB) stars are more concentrated in 
the central region than old RGB stars. 

A new, independent data set of metallicities for Leo I stars, 
also based on the CaT method, was obtained by us at the ESO 
VLT using high signal-to-noise FORS2 spectra. The new sam- 
ple has negligible overlap with the previous data sets, thus ef- 
fectively increasing the number of Leo I stars with direct metal- 
licity measurements. Using the new data, this paper provides an 
independent determination of the MDF of Leo I and further con- 
straints on its evolution, based on an analysis of metallicity and 
age gradients and the age-metallicity relation. 



2. Observations and reduction 

2.1. Target selection 

Our targets were selected from the colour-magnitude diagram 
(CMD) of Leo I. For the central re gion of the galaxy , we relied 
upon B, V-band photometry from IHeld et ail (|2000|) . obtained 
with the EMMI instrument at the NTT telescope at the ESO La 
Silla Observatory. For the stars in the outer regions, we used 
the BV photometry ori ginally obtained for a study of RR Lyrae 
variable stars in Leo I (He ld et al. I l200ih . based on observations 
carried out with the Wide Field Imager at the 2.2m ESO-MPI 
telescope. 

The spectroscopic targets were selected among the bright- 
est RGB stars of Leo I, down to 1 mag below the RGB tip. We 
avoided any colour constraints that might bias the age/metallicity 
distribution. After mask design (for which we were guided only 
by geometric constraints) we were left with 61 targets in 4 
masks. 

The identifiers, coordinates, and BV photometry of the stars 
in our final sample (excluding a few targets with too low S/N 
ratio to allow any measurements) are listed in Table [1] and the 
targets are shown in Fig.Q] The location of the target stars in the 
CMD of Leo I is shown in Fig. [2] 



2.2. Observations 

The observations were carried out in service mode in two runs 
between May 2002 and December 2003 using FORS2, the 
multi-mode optical instrument mounted on the Cassegrain fo- 
cus of the Yepun (VLT-UT4) 8.2m telescope at the ESO Paranal 
Observatory. We used FORS2 in MXU mode with the 1028z+29 
grism and the OG590+32 order-blocking filter. With this setup, 
the spectral coverage is approximately 7700 A to 9500 A, with 
a dispersion 0.85 A pixel -1 . The selected targets were observed 
with 4 masks using 0780 slits. For each Leol mask, two spec- 
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Fig. 2. Target stars in the colour-magnitude diagram of Leo I. 
The filled squares represent red giant stars in Leo I, the diamonds 
are the 4 metal-poor red giant members, while crosses mark non- 
members. 

Table 2. Log of the observations. 



Field 


Night 


tap (sec) 


NGC 4590 


6 May 2002 


15 + 60 


NGC 5927 


6 May 2002 


60 + 300 


NGC 6171 


6 May 2002 


10 + 60 


NCG 5904 


3 May 2003 


60 + 300 


NGC 6397 


5 May 2003 


60 + 300 


NGC 6528 


5 May 2003 


60 + 300 


NGC 4372 


24 May 2003 


60 + 300 


NGC 6752 


24 May 2003 


60 + 300 


Leo I field U 


4 May 2003 


2 x 2870 


Leo I field R 


20,21 Dec 2003 


2 x 2870 


Leo I field D 


22 May 2003 


2 x 2870 


Leo I field L 


16,25 Jun 2003 


2 x 2870 



tra were taken. The observing log and exposure times are given 
in Table [2] We observed also RGB stars in 8 Galactic globular 
clusters (GCs) in a wide range of metallicity with the same in- 
strumental setup in order to calibrate our measurements onto a 
known metallicity scale. Two different exposure times were used 
(long and short exposure) to prevent saturation of the brightest 
RGB stars. The ph otometry of the GC star s was taken from the 
data compilation of lRutledge et al.l (Il997bl) . 

2.3. Data Reduction 

The basic reduction of multi-object spectra was performed us- 
ing standard procedures in IRAF Q Bias and flat-field correc- 
tions were applied to all images using the ccdproc task. Due 

1 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under cooperative agreement with the 
National Science Foundation. 




8500 8550 8600 8650 8700 

x(A) 

Fig. 3. Examples of normalised, background-subtracted spectra 
of RGB stars in Leo I (stars #17, #19, #54, from top to bottom). 
The spectra are shown on a rest-frame wavelength scale and ver- 
tically shifted for clarity. 

to the long exposure times, the scientific frames for Leo I con- 
tain a large number of cosmic ray hits. These were effectiv ely 
cleaned using the IRAF program lacos (Ivan Dokkumll200l on 
the bias-subtracted, flat-fielded images. The multi-object spec- 
tra were extracted with the apall task and wavelength calibrated 
using HeNeAr lamp exposures taken at the end of each night. 

The two spectra taken for each Leo I target were combined 
to increase the S/N ratio. However, we also retained the two in- 
dividual spectra to estimate the uncertainties in the wavelength 
calibration and line strengths. Finally, the continuum was nor- 
malised in the region between 8400 A and 8800 A, by exclud- 
ing in the process the Ca n and other relatively strong absorption 
lines. Typical sky-subtracted spectra of three Leo I stars with dif- 
ferent metallicities are shown in Fig. [3] 

The average signal-to-noise ratio per pixel was calculated 
from the rms of the combined spectra in two wavelength win- 
dows free from strong spectral features, 8580-8620 A and 87 10- 
8750 A. The values were checked against those measured on the 
raw spectra, and found to be consistent within a few percent. The 
S/N ratio, listed in Table Q] is > 20 for all stars in our sample, 
with a mean value of ~ 50. 

3. Radial velocities and membership 

Radial velocities were measured for target stars to establish their 
membership. Only the two reddest lines were used to this pur- 
pose, since the bluest line of the Ca triplet is weak and, for the 
systemic velocity of Leo I, overlapping with a strong sky line. 
The line wavelengths were obtained from the central values A m 
of the fitted profile (see Eq. Q] below), and compared with the 
laboratory air wavelengths 8542.09 A and 8662.14 A. A radial 
velocity was measured for each individual spectrum of each star 
by combining the measurements of the two lines /I8542 and /18662- 
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Fig. 4. Heliocentric radial velocities of stars in the Leo I field. 
The distribution is fitted with a Gaussian centred at 271 km 
s _1 and with dispersion 13.7 km s _1 . The vertical dashed lines 
represents the 3<x limits used to select members of Leo I. 
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Fig. 5. Upper panel: correlation between two independent EW 
measurements on individual spectra of Leo I stars. The dashed 
lines show the 3<x interval. Lower panel: histogram of the dif- 
ferences between the two measurements, fitted by a Gaussian 
profile with cr^ w = 0.44 A. 



Then, the radial velocity was calculated as the mean of the two 
values independently measured from the single spectra. The re- 
sults, corrected to heliocentric velocities using the rvcorr task, 
are given in Table Q] The distribution of radial velocities, shown 
in Fig. [4] is well fitted by a Gaussian function centred at 271 km 
s _I with a dispersion 13.7 km s . All but three stars have he- 
liocentric radial velocities within 3<x of the peak, and therefore 
are considered members of Leo I (members have ID < 100 in 
Tabled). 

Radial velocity errors (given in Col. 7 of Table Q]) were es- 
timated by calculating the differences between radial velocities 
measured on the two individual spectra. The distribution of the 
differences is approximately Gaussian with a mean -6.4 km 
s _1 and dispersion 10.4 km s _1 . The mean velocity difference 
provides a good estimate of the mask-to-mask systematic er- 
rors, which include wavelength calibration errors and the cen- 
tring errors of the targets on the slitlets, while the standard de- 
viation represents a combination of random errors. We assume 
cr Y = cr Vl _ V2 / y2 = 7.4 km s _1 as a good approximation to 
the standard error of the radial velocity measured on the com- 
bined spectrum. Therefore our systemic radial velocity of Leo I 
is 271 ± 6.4 (systematic) ± 7.4 (random) km s _1 . Our error es- 
timates are consistent with the 0.85 A pixeL 1 spectral resolu- 
tion provided by our instrumental setup, which corresponds to 
a resolution in radial v elocity of ~ 30 k m s^'pixeF 1 in the 
CaT wavelength range. iKoch et a i1 (l2007hl) measured 284.2 km 
s _I with a velocity dispersion of 9.9 km s" 1 , while 2 82.6 ± 9.8 
km s~ 1 was the value meas ured by iBosler et alJ (|2007 ) . Most re- 
cently, [MatecTeTal] d2008l) obtained a mean heliocentric veloc- 
ity 282.9 ± 0.5 km s and a dispersion 9.2 ± 0.4 km s _1 from 
echelle spectroscopy of 328 Leo I members. The systematic dif- 
ference of ~ 10 km s _1 between our mean velocity estimate and 
previous results is consistent with pointing errors. In the CaT 



wavelength range, a velocity of 10 km s 1 corresponds to a shift 
of ~ 1/3 pixel on the detector (see above), or U.'OS on the sky. 



4. Equivalent widths and metallicity 

4. 1 . Equivalent width measurements 

We measured the equivalent widths (EWs) of CaT lines in the 
spectra of target stars in Leo I and the calibrating GCs as fol- 
lows. We first normalised the spectra over a wavelength in- 
terval encompassing, fo r each line, the side bands defined by 
lArmandroff & Da Costal d 199 lb - Equivalent widths were then 
measured for the two stronger CaT lines in the co-added spec- 
tra by fitting a model profile over the lin e central bandpass es as 
defined by the same authors. Following ICole et"aT] ( |2004|) . the 
fitted model is the sum of a Gaussian and a Lorentzian profiles 
with a common line centre /L, 



F(A) = 1 - A G exp 



(A-A m ) 2 ' 




2o- 2 





r 1 



+ 1 



(i) 



with the best-fit parameters determined using a Levenberg- 
Marquardt least-squares algorithm (coded in an idl procedure 
by C. Markwardu). The Can line strength was then defined as 
the unweighted sum of the two equivalent widths, 



ZW = £W 85 42 + £W 86 62 



(2) 



The sum equivalent widths are given in Table |3]for the red giants 
in Leo I and in Table|4]for the RGB stars in the template GCs. 



http://cow.physics.wisc.edu/~craigm/idl/idl.html 
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Table 4. Observed stars in template globular clusters. 



Fig. 6. The sum of the equivalent widths of the two most reliable 
CaT lines plotted as a function of the magnitude difference from 
the HB level. The calibration globular clusters are represented 
by open symbols (coded with different colours in the electronic 
edition of the journal). The metallicity of the clusters increases 
from the bottom to the top (see Table 13). The straight lines are 
best fits to the EWs for each cluster, assuming a common best-fit 
slope. Filled squares refer to Leo I RGB stars. 

To estimate the equivalent width measurement errors, we 
also measured the CaT line strengths independently on the indi- 
vidual spectra of each star in Leo I. A comparison of the two YW 
measurements is shown in Fig. |5j where they appear to be well 
correlated. The histogram of the differences A^jy = (LWo - SWi) 
(shown in Fig.[5j lower panel) is well fitted by a Gaussian func- 
tion centred at (A^w) = 0.03 and with a standard deviation 
o"aw = 0.44 A. Since the two individual spectra have compara- 
ble S /N ratio, we adopt cr-^w = cr^wl V2 = 0.31 A as our error 
estimate for YW measured on the combined spectrum. For com- 
parison, the half-range ej.w - \Azw/2\ is listed for each Leol 
star in Table [3] 

4.2. Reduced EW 



According to the CaT method (lArmandroff & Da Costal fl99ll) . 
the gravity and dependence of CaT lines is accounted for 
by introducing a linear correction to the line strengths which de- 
pends on the star's luminosity, that is 



W = YW + B(V - Vhb) 



(3) 



where B is a constant and (V - Vhb) is the difference be- 
tween the magnitude of the star and the horizontal branch 
(HB) in the V band. In globular clusters, this reduced equiv- 
alent width W was found to be well correlated with metal- 
licity (lArmandroff & Zinnll 19881: lArmandroff & Da Costalll991t 
iRutledge et alj[l997bl) . This provides the empirical basis for the 
validity of the CaT method. Applied to composite stellar popu- 
lations, the method is less straightforward and has been widely 



Cluster 


ID 


V 


B-V 


YW(A) 


M5 


11-45 


14.75 


0.82 


4.62 




11-50 


13.92 


0.96 


5.12 




11-51 


14.05 


0.96 


4.90 




11-80 


14.31 


0.91 


4.89 




11-74 


13.82 


1.01 


5.05 




1-2 


13.87 


1.02 


5.24 




1-50 


13.91 


0.97 


5.06 




1-61 


13.37 


1.17 


5.40 




1-68 


12.37 


1.52 


6.25 




1-71 


13.01 


1.29 


5.57 


NGC 4372 


20 


12.88 


1.57 


3.78 




14 


14.29 


1.29 


2.88 




13 


12.72 


1.73 


4.13 




10 


13.82 


1.29 


3.00 




95 


14.48 


1.30 


2.72 




91 


14.45 


1.33 


2.81 




89 


14.49 


1.29 


2.83 




141 


12.93 


1.65 


4.01 




74 


14.17 


1.40 


3.09 




76 


14.18 


1.30 


3.16 




77 


14.19 


1.28 


3.07 


NGC6171 


62 


13.97 


1.62 


5.83 




100 


14.21 


1.40 


5.50 




I 


13.89 


1.46 


5.66 




F 


13.39 


1.70 


6.21 


NGC 6397 


328 


12.07 


0.93 


3.19 




326 


12.78 


0.89 


2.89 




337 


12.58 


0.90 


2.83 




343 


11.42 


1.13 


3.34 




361 


11.67 


1.08 


3.18 


NGC 6528 


R2-8 


15.79 


1.89 


6.80 




Rl-42 


16.46 


1.62 


6.34 




R2-41 


16.30 


1.64 


6.44 


NGC 6752 


4 


13.70 


0.85 


3.86 




8 


11.96 


1.05 


4.99 




28 


13.17 


0.89 


4.18 




29 


11.79 


1.17 


4.98 




30 


12.15 


1.12 


4.83 


NGC 5927 


133 


14.75 


1.97 


6.45 




372 


14.66 


2.11 


6.28 




335 


14.44 


1.94 


6.63 




190 


14.29 


2.02 


6.71 




65 


14.64 


1.92 


6.58 


NGC 4590 


144 


12.80 


1.29 


3.72 




239 


14.19 


0.87 


2.88 




1172 


15.03 


0.85 


2.60 




30 


14.15 


0.87 


2.64 




74 


14.59 


0.84 


2.36 




119 


13.62 


0.95 


2.90 



Notes. The IDs, magnitudes, and colo urs of stars are t hose gi ven in the 
original photometry papers quoted by Rutledge et al. ( 1997a). 



discussed in the recent literature. We will return to this point later 
on. 

In Fig. [6] we plot the sum of equivalent widths YW versus 
V - Vhb for all stars with good S/N spectra, both in Leol and 
the template GCs. The magnitude of the HB (of old stars) in 



M. Gullieuszik et al.: New constraints on the chemical evolution of the dwarf spheroidal galaxy Leo I from VLT spectroscopy 



Table 5. Parameters for the calibration Galactic globular clus- 
ters. 



Cluster 


[Fe/H] Z w 


[Fe/H]cG 


[M/H] 


Vhb 


NGC 6528 


-0.23 


-0.10 


-0.03 


17.10 


NGC 5927 


-0.31 


-0.46 


-0.37 


16.60 


NGC 6171 (M107) 


-0.99 


-0.87 


-0.70 


15.70 


NGC 5904 (M 5) 


-1.40 


-1.11 


-0.90 


15.06 


NGC 6752 


-1.54 


-1.42 


-1.21 


13.70 


NGC 6397 


-1.91 


-1.82 


-1.65 


12.87 


NGC 4372 


-2.08 


-1.94 


-1.74 


15.30 


NGC 4590 (M68) 


-2.09 


-1.99 


-1.81 


15.68 



Table 6. Metallicity calibrations. 



Type 



Calibration 



quadratic 


[Fe/H] ziv = 


0.088 W' 2 


-0.184 W 


-2.079 




[Fe/H] CG = 


0.072 W' 2 


-0.076 W 


-2.122 




[M/H] = 


0.051 W' 2 


+0.056 W 


-2.125 


linear 


[Fe/H] ZW = 




0.359 W 


-2.845 




[Fe/H] CG = 




0.391 W 


-2.806 




[M/H] = 




0.395 W 


-2.628 



Leo I, V HB = 22.60, is from lHeld etai1d2001l) while fo r the tem- 
plate GCs, Vhb was taken from Rutledge et alj dl997al) (listed in 
Table [5] together with clusters' metallicities). For all the globu- 
lar clusters, our CaT line strengths define clean, well separated 
linear sequences generally consistent with a constant slope and 
having different, metallicity-dependent zero points. In this di- 
agram, the Leo I stars show quite a large dispersion, although 
most of them are located between the sequences of NGC 6397 
and M5. This spread in CaT line strengths is real, being larger 
than the typical measurement error cr-^w ~ 0.3 A. By assuming 
a common slope and lumping together the data for all globular 
clusters, we derived a slope f3 = 0.627 ± 0.021. This value can 
only be compared with previou s results that use the same defini- 
tion of YW. This is the case for lTolstov et al.l d200ll) . who found 
ft = 0.64 + 0.02, in agreement with our result. 

4.3. Metallicity calibration 

Using our reduced CaT equivalent widths W and the published 
metallicities for the Galactic globular clusters, we re-determined 
the calibration relations between W and m etallicity on 3 differ- 
ent abundance scales: the [Fe /H] scales of IZinn & West! (1 19841) 
and lCarretta & Grattonl(ll997l) and the global metallicity [M/H], 
as defined bv lSalaris et a l. ( 1993). The metallic ities of the GCs 
(Table [5]l were taken from iFerraro et al.l (1 19991) . except for the 
m etal-rich cluster NGC 6528, for which the more recent results 
of IZoccali et ail d2004l) were adopted ([Fe/H] = -0.1, [ar/Fe] = 
0.1). The global metallicity of N GC 6528 ([M/H] = - 0.03) was 
calculated using the relation from lSalaris et alj (1 1 993b : 



[M/H] =[Fe/H] + log(0.683/„ + 0.362) 
log f a = [ar/Fe] 



(4) 



Figure [7] shows the W -metallicity relations for the three 
scales, along with quadratic fits to the whole dataset and lin- 
ear fits to the metal-poor and intermediate globular clusters. 
The quadratic relations provide a better fit to the GC metal- 
licities over the whole metallicity range of template GCs. 
The curvature is driven by the data for two most metal- 
rich globular clusters, consistently with a fall in sensitiv- 
ity of the Can index at high metallicity. Previous studies 




Fig. 7. The metallicity of the reference globular clusters against 
the reduced EW of CaT lines, for three adopted metallicity 
scales. The solid lines are quadratic fits, while the dashed 
straight lines are linear fits obtained for clusters with [Fe/H]< 
-0.6. 



which included metal-rich GC s also found quadratic relations 
dArmandroff & Da Costal [1991: iDa Costa & Amandrofllll995t 
ICarretta et all 1200 U iBosler et aill2007l) . Line ar relations have 
been proposed by other studies (most re cently, ICole et al]|2004l : 
iKoch et al] l2007bt ICarrera et alj l2007t) using metal-rich open 
clusters to constrain the metal-rich end of the W - [Fe/H] re- 
lation. A full discussion of the behaviour of CaT line strengths 
against metallicity is beyond the scope of this paper, and will be 
presented in a future paper along with a large dataset of calibrat- 
ing globular clusters. For our data, a linear relation indeed pro- 
vides a good fit for stars less metal-rich than the template cluster 
NGC 6171 ([M/H] = -0.70) (Fig. [7). In the case of a metal-poor 
system such as Leo I, the linear and quadratic relations give sim- 
ilar results except for the most metal-poor stars. Our calibration 
is presently quite uncertain near the metal-poor end, being based 
on one globular cluster (NGC 4590). For this cluster. iPritzl et al.l 

give a lower metallicity ([Fe/H] 2.3) than that adopted 

in Table [5] yielding a better agreement with our linear calibra- 
tion. 

The main source of error on [Fe/H] (or [M/H]) is the uncer- 
tainty on the measured equivalent width YW, since other sources 
of error, such as photometric errors for stars on the upper RGB 
of Leo I, the error on the HB level, or even the uncertainties 
associated to the fit parameters of the calibration relations, are 
negligible compared to the YW measurement errors. A metallic- 
ity uncertainty can be computed for each star by error propaga- 
tion using the values of e zw in Table [3] and the calibrations in 
Table [6] However, a more meaningful metallicity uncertainty is 
obtained using cr w , as our estimate of the measurement error. For 
the quadratic [M/H] calibration in Table|6j crw = 0.3 1 A implies 
a [M/H] error cr elT 0.14 dex for stars with [M/H] = -1.20. 

The sources of uncertainty related to the CaT method itself 
are more difficult to quantify and predict. All traditional calibra- 
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tions refer to Galactic globular clusters, which are simple and 
nearly coeval old stellar populations, and the applicability of 
these calibrations to complex stellar populations is not obvious. 
In our case, Leo I stars are on average several Gyr younger than 
those in GCs. At a given luminosity, a star in Leo I has a differ- 
ent mass from a GC star with the same metallicity and V - Vhb- 
However, recent studies have shown that the CaT method can be 
used for com plex stellar popul a tions younger that those in glob- 
ular c l usters dCole et al J 120041: iPont et al.ll2004t IBattaglia etail 



Table 7. Mean metallicity and standard deviation of red giants 
in Leo I. 



120081) . IBattaglia et al.l d2008l) have compared a linear metallic- 
ity calibration that uses the CaT of RGB stars in two dSph 
(Sculptor and Fornax) with spectroscopic [Fe/H] values obtained 
from their high-resolution studies. The metallicities are in good 
agreement, although with some residual trends of about 0.1-0.2 
dex, in the range -2.5 < [Fe/H] < -0.8. These studies conclude 
that for ages older than 2.5 Gyr, the CaT line strengths are little 
affected by age, and suggest that the overall uncertainty related 
to age effects is < 0.2 dex. 

4.4. A new metallicity scale? 

Most CaT metallicity measurements in nearby galaxies use 
a [Fe/H] sc ale based on observations of Galactic globular 
cluste rs (e.g. IPont et al.ll2004t ICole et all 120041: IBattaglia et ail 
120061: iKoch et al.ll2006l) . However, the relative abundances of a- 
elements (including Ca) in Local Group dwarfs are on average 
lower than in the Milky Way halo stars and GCs (|Shetrone et al 



2001 



120031: iTolstov et al.ll20Q3b iGeisler et alJl2005t iPritzl et al 



2005). To overcome this problem. iBosler et al.l (12007) proposed 
a new calibration of CaT lines against the [Ca/H] abundance, 
based on high-resolution spectroscopy of Galactic star clusters. 
In their hypothesis, the [Ca/H] calibration is less affected by the 
difference in [Ca/Fe] abundance ratios between red giant stars in 
globular clusters and dwarf spheroidal galaxies. However, the 
strength of CaT lines is also determined by other parameters 
(gravity and T e ff) in addition to Ca abundance. A comparison of 
Ca abundances derived from CaT lines with the results of high- 
res olution abundance me asurements for stars in two dSph galax- 
ies ( IBattaglia et al J 120081) shows that, while the CaT lines trace 
both Ca and Fe, their dependence on Fe abundance is stronger. 
Similarly, [Ca/H] ratios deri ved from CaT lines fo r stars in Leo II 
dSph and globular clusters dShetrone et al.ll2009l) systematically 
differ from those obtained from mid-resolution synthetic spectra, 
with a residual trend that is a function of metallicity. Since the 
effective temperature of red giant s in globular cluste rs is driven 
by their global metallicity [M/H] ( ISalaris et al.ll 19931) . an empir- 
ical metallicity ranking based on CaT and a global metallicity 
[M/H] scale, proposed here for the first time, appears to be the 
most empirically sound. In fact, [M/H] (or, equivalently, Z) takes 
into account the abundances of both the ff-elements and Fe. Our 
data are therefore calibrated using the [M/H] calibration in addi- 
tion to the common [Fe/H] scales. 



4.5. A concluding remark 

We conclude this section with a consideration that ought to be 
kept in mind throughout all the following discussion. While 
we give the metallicity of the stars in three flavours ([Fe/H] zw , 
[Fe/H] CG , and [M/H]), this does not imply that we are determin- 
ing the three metallicity parameters at the same time. The only 
observable quantity is the reduced equivalent width W. The cali- 
bration of W in terms of metallicity relies on the assumption that 
W is correlated with metallicity, i.e. a star in Leo I has the same 



Scale 
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mean 


cr 


[Fe/Hf w 


linear 


-1.53 


0.17 




quadratic 


-1.55 


0.21 


[Fe/H] CG 


linear 


-1.37 


0.18 




quadratic 


-1.41 


0.21 


[M/H] 


linear 


-1.18 


0.19 




quadratic 


-1.22 


0.20 



metallicity as a star in a GC with the same W. The key questions 
are: what are the real drivers that determine the strength of the 
CaT lines? Do two stars with the same iron-to-hydrogen ratio 
but different a-elements composition ha ve the same W' l Som e 
of these effects have been discussed by Batta glia et al.l d2008l) . 
and are further addressed by a large observational program by 
our group whose results will be presented in future papers. 



5. The metallicity of Leo I stars 

5. 1 . The observed metallicity distribution 

In Fig. |8] we show the metallicity distribution of Leo I red giant 
stars as de rived from our data using both the [Fe/H] metallic- 
ity scale of ICarretta & Grattonl dl997t) and the [M/H] scale. The 
parameters of the distribution (mean and standard deviation, ex- 
cluding the 4 most metal-poor stars in Fig. [8 1 are given in T able[7] 
for both scales, along with the results on the lZinn & Westl dl984l) 
scale for ease of comparison with previous literature. The metal- 
licities of individual Leo I stars are listed in Table [3] Using the 
CG97 scale, the distribution is centred at [Fe/H] ^ -1.4 with a 
standard deviation cr^/u] - 0.2. For the [M/H] calibration, the 
average is [M/H] ^ - 1 .2 with the same scatter. 

The results obtained from the linear and quadratic calibra- 
tions are very similar in all cases, as expected since most of the 
Leo I stars have metallicities lower than [Fe/H] = -1.0. The 
choice of the linear or quadratic relation only affects the metal- 
licity of 4 metal-poor stars, having [M/H]< -2. In this range, the 
calibration is extrapolated beyond the most metal-poor globular 
cluster, which makes the metallicity of the 4 stars quite uncer- 
tain and dependent on the adopted calibration. In the case of the 
CG97 calibration, which yields the lowest extrapolated values, 
the 4 stars have -2.6 < [Fe/H] < -2.2. Visual inspection of the 
targets on a Leo I image indicates normal star-like profiles (i.e. 
no blends). We therefore conclude that, while a few stars may 
have low metallicity, there is so far no evidence of extremely 
metal-poor stars i n Leo I. Spectral synthesis methods (see, e.g., 
iKirbv et al.ll2008l) will be used in a future paper to obtain more 
secure metallicity estimates for these metal-poor stars from a 
different spectral interval. 

Our determinat ion of the metallic ity of Leo I agrees well 
with the results of IBosler et al.l (120071) ([Fe/H] = -1.34) and 
IKoch et all d2007bl) ([Fe/H] = -1.31), in particular when a lin- 
ear calibration is used as in the previous papers. On the other 
hand, our [M/H] value s are in better agreement with the [Ca/H] 
results of IBosler et al.l d2007l) . Using 9 stars in co mmon with 
IBosler eta l. (2007) and 5 stars in common with IKoch et ail 
(I2007bl) . we compared the metallicities star-by-star. The mean 
differences are A[Fe/H] CG = -0.04 + 0.14 (rms of the sample) 
and A[Fe/H] CG = -0.17 ± 0.11, respectively, in agreement with 
the shifts between the mean values of the MDFs. The rms values 
are consistent with our measurement error (see next section). 
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Fig. 8. Metallicity distribution of Leo I stars using a linear (left panels) or quadratic (right panels) calibration re lation. In the upper 
panels we plot the MDF obtained with the [M/H] calibration. The metallic ity distributions on the [ Fe/H] scale of ICarretta & Grattonl 
dl997l) are shown in the lower panels, together with previous results from lBosler et al.l d2007l) and lKoch et al.l d2007bl) . 



5.2. Intrinsic metallicity dispersion and clues on the evolution 
of Leo I 

The observed MDF in Fig. [8] is the convolution of the intrin- 
sic metallicity distribution of stars in Leo I and the measure- 
ment errors. The real abundance spread can be estimated by 
adopting a Gaussian model for the intrinsic MDF, which yields 
°"obs ~ <T o + °eir' where <x 0BS is the observed metallicity disper- 
sion, <xo is the intrinsic dispersion, and cr en - is the measurement 
error. 

If we adopt a quadratic [M/H] calibration and a typical mea- 
surement error cr elT 0.14 dex, the measurement scatter largely 
contributes to the observed metallicity dispersion. The observed 
scatter implies an intrinsic metallicity dispersion cr[ Fe / H ] o = 0.14 
for the CG97 scale and ct[m/h],o = 0.08 for [M/H]. The intrinsic 
abundance dispersion of Leo I stars is therefore very small, even 
smaller than previously thought, and this happens in spite of the 
relatively wide range of ages of the stellar populations. This is 



an important constraint to the chemical evolution across the life 
of the galaxy. 

In order to model the metallicity distribution and the chemi- 
cal evolutio n of Leo I, detailed models have to be put forth (such 
as those of lLanfranchi & Matteuccill2007l) properly taking into 
account the chemical and dynamical evolution of the galaxy. 
However, some order-of-magnitude physical insight on the evo- 
lution of Leo lean already be obtained using basic considera- 
tions. The metallicity distribution of RGB stars in Leo I is com- 
pared in Fig. [9] with a simple closed-box model with a l ow ef- 
fectiv e yield consistent with a continuous loss of gas (e.g. JPagel 
119971) . In order to reproduce the peak of the observed MDF, we 
have to adopt an effective yield y = O.O25Z , and y = 0.040Z G , 
for the distributions based on the [Fe/H] and [M/H] metallicity 
scales, respectively. This is clearly much lower than the valu e 
found in the solar vicinity (y = 1.2Z Q ; e.g. lPortinari et al]|2004l) . 
in a way consistent with the loss of metals driven by a galactic 
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Fig. 9. The MDF of Leo I RGB stars on the [M/H] metallicity 
scale, compared with a simple model with a low effective yield 
(dashed line) and a model with a prompt initial enrichment (solid 
line). 

wind dHartwickl 1 976t iPagell 1 99% . Still, even allowing for a gas 
outflow, the number of metal-poor stars is largely overestimated 
by the simple model, as shown in Fig. [9] The fit is considerably 
improved by assuming a prompt early enrichment with an ini- 
tial metallicity Zq = 0.02 Z Q ([Fe/H] = —1.7) (continuous line 
in Fig. [9]). Although very simplistic, this conclusion agrees with 
the finding that the metal- poor tail of the MD F in 4 Local Group 
dwarf spheroidal galaxies dHelmi et al.l2006l) is significantly dif- 
ferent from that of the Galactic halo, lacking stars below [Fe/H] 
= -3. What this "toy model" tells us is that the narrow MDF of 
the Leo I stars can be understood as a combination of fast enrich- 
ment from an initial generation of stars, and subsequent loss of 
metals through outflows. This situation is common among Local 
Group d warfs, but the MDF of Leo I is the narrowest observed to 
date (c f. iTolstov et aUbOQlHPont et aDl2004t iKoch et alJl2006L 
l2007ah . 

5.3. Radial metallicity gradients 

Radial variations in the stellar populations are common in 
the dwarf spheroidals of the Local Group, where the younger 
and more metal-rich populations are often concentrated to- 
ward the galaxy centre dHarbeck et al.l200ll;ISaviane et alj200lt 
iPont et alj|2004t ITolstov et al.ll2004tlKoch et al.ll2006l) . Leo I re- 
mains one of the few dSph's showing little eviden ce of a popula- 
tion gradient dHeld et alJ l2000l: lKoch et"aT]l2007bl). In particular, 
the spectroscopic investigation of lKoch et al. I d2007bl) . extending 
to quite large radial distances, did not detect a significant metal- 
licity gradient. 

Our new spectroscopic sample of red giants allows us to fur- 
ther search for radial variations in the metallicity of Leo I stars. 
As the radial coordinate, we have adopted the semi-major axis r 
of ellipses passing through the projected sky position of each 
star. The ellipses have the centre at 10*08 m 2Kl, +12°18'23" 
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Fig. 10. Metallicities of Leo I stars on the [Fe/H] scale of CG97, 
plotted against the elliptical radius (s ee text). Filled dot s: data 
in this pa per; open squares: d ata from lBosler et al.l d2007l) ; open 
triangles: IKoch et all (l2007bb : circles with error bars: our data, 
binned in 1 '5 bins. The error bars of the binned data represent 
the abundance scatter (lcr) in each bin. The crosses are the 4 
metal-poor stars in our sample. The typical errors of each study 
are shown in the upper right corner of the plot. The solid line is 
a fit to our (binned) data, while the dashed line represents a fit to 
all available spectroscopic data. 



(J2000) a nd a fixed position angle and ellipticity (PA = 79°, 
e = 0.21; llrwin & Hatzidimitrioull 19951) . 

In Fig. [T0l the metallicities of Leo I stars in our sample (see 
Table [3]) are plotted against the distance from the centre. To di- 
rectly compare our data with results from previous studies, we 
used the metallicity obtained from the cali bration in terms of 
[Fe/H] CG . Data from lBosler etafl d2007l) and lKoch et all d2007bl) 
were shifted to account for the small differences in mean metal- 
licity (of the order 0.1 dex or less) between the MDF's (Fig. [8j. 
A linear relation was fitted both to our data alone (solid line) 
and to all available metallicity measurements in the literature 
(dashed line), excluding stars with [Fe/H] < -2. 

The fit to our new FORS2 data yields a radial gradient of 
-0.02 dex arcmin -1 , or -0.27 dex Kpc -1 . In our sample, stars 
more metal -rich than [Fe/H] = -1 .3 are only found in the central 
region of Leo I, with a < 5'.5. In contrast, the fit to the merged 
spectroscopic sample suggests a r adially cons t ant me tallicity, in 
agreement with the conclusions of lKoch et al.l(l2007bl) . To quan- 
tify the gradient, we have used a Kolmogorov-Smirnov test to 
compare the metallicity distributions of stars with V.2 < a < 5' 
(inner sample) and 5' < a < 8'2 (outer sample), considering 
only the radial interval covered by our data. The two metallic- 
ity distributions are similar in shape, with the MDF in the inner 
region peaked at higher metallicity. The hypothesis that the in- 
ner and outer sample are drawn from the same parent population 
can be rejected at a 90% level using our data, and only at a non- 
significant 62% level using all spectroscopic data. 

We conclude that, while our data provide a hint of a weak 
radial metallicity gradient in Leo I, the statistical significance of 
this result is at present low. More stars need to be observed, par- 
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Fig. 11. The age-metallicity relation of Leo I RGB stars in our 
sample, on the [M/H] scale. The error bars in age represent 
the first and third quartile of the confidence intervals obtained 
through Monte Carlo realisations (see text for details). For the 
metallicity errors the representative value discussed in Sect. 14.31 
is adopted. Also shown are the metallicities from high-resolution 
spectroscopy of two Leo I RGB stars from lTolstov et al] d2003l) 
(big filled circles). The side histograms are the marginal distri- 
butions in metallicity and age. The sol i d line in the top panel 
represents the SFH derived by iDolphinl (|2002|) from HST pho- 
tometry, normalised to the total number of stars in our sample. 

ticularly in the outer region of the galaxy, before definite conclu- 
sions can be drawn. 

6. The age of Leo I stars 

6.1. The age-metallicity relation 

With the stellar metallicities of Leo I stars known from spec- 
troscopy, ages could be estimated by comparing stars' loca- 
tions in th e CMP with a grid o f theoretical isochrones (the 
models of iPietrinferni et al.|[2004 were used to this purpose). 
Absolute magnitudes and dereddened colours were computed 
adopting a colour excess Eb-v - 0.04, a total-to-differential ex- 
tinction ratio Ry = Ay/Es-y - 3.1, and a true distance modulus 
(m - M) = 22.04 dHeld et alj|2009b . 

The stellar ages were interpolated in two steps. First, we 
used a set of theoretical isochrones of fixed age and different 
metallicities to find, for a star of given age and luminosity, a 
metallicity-colour relation and (from the known colour) an in- 
terpolated metallicity. This step was repeated for all model ages, 
yielding for each data point in the CMD a set of theoretical 
age-metallicity pairs each consistent with the star's magnitude 
and colour. This age-metallicity look-up table (spanning the full 
range from old, metal-poor stars to young, metal-rich stars) al- 
lowed us to compute an interpolated age for each star from its 
spectroscopic metallicity. We chose the [M/H] scale for the input 
value, as the most directly related to the mass fraction of metals 
(Z) used in stellar models. For a number of stars, ages could not 



be derived because the observed colour and/or magnitude were 
outside the range covered by the isochrones. 

The method was checked against Galactic globular clusters 
with ages given by the literature. In particular, a small correc- 
tion was applied to the isochrone colours so as t o yield a correct 
age fo r NGC 5904 (M 5), assumed to be 12 Gyr dSandquist et al.l 
1996) and the closest in metallicity to Leo I among the clusters 
listed in Table [5] Thus, our ages for Leo I stars are essentially 
referred to M5, which is in our view the most correct approach 
given the considerable uncertainties in the isochrone colours. 

The resulting ages are listed in Table [3] along with their un- 
certainties, estimated as follows. For each star, we performed a 
set of 100 Monte Carlo experiments by randomly varying the in- 
put quantities in intervals consistent with their uncertainties. We 
adopted a standard error 0.02 mag in the V magnitude, 0.05 mag 
in (B - V) colour, and 0.15 dex in [M/H]. The latter was chosen 
conservatively large to account for the inherent uncertainties in 
the metallicity scale. The median and quartiles of the age dis- 
tributions of randomly generated "stars" corresponding to each 
observed stars are listed in Table [3] In general, the uncertainties 
are of the order 50%, which reasonably reflects the large uncer- 
tainties in the process. 

The age-metallicity relation derived from our data is shown 
in Fig. QjJ The chemical evolution of Leo I seems to be very 
slow, in accord with the narrowness of the MDF If the 4 
metal-poor stars are excluded, there is a trend for stars younger 
than 5 Gyr to be on average more metal-rich by about 0.2- 
0.3 dex . Similar conclusion s were drawn, from a different data 
set, by iBosler et al.l d2004l) . The scatter in the age-metallicity 
relat ion appears to be smaller than observed in other galax- 
ies dBattaglia et alj l2006t iTolstov et al.l I2003L and references 
therein). Our data are in agreemen t with the results fr om high- 
resolution spectroscopy for 2 stars dTolstov et alj|2003l) . 

We can use our age determinations also to obtain a SFH of 
Leo I. This can be done since our target selection was designed to 
avoid any bias in age and/or metallicity (see Sect. [2]); the age dis- 
tribution of our target stars is then proportional to the SFH. Our 
age measures are in agreement with the SFHs derived by HST 
photometry. In the u pper panel of F ig. [TTIour age distribution 
is compared with the IDolphinl d2002l) SFH, showing only small 
differences that can be explained by statistical fluctuations. 

6.2. Radial distribution of stellar ages 

To complete our analysis of the population gradients, we inves- 
tigated the possible presence of a radial variation in the age of 
Leo I stars. Figure[T2]shows the ages plotted against the elliptical 
distance. While old stars are found at all radii, young stars appear 
to be concentrated at small distances from the centre. To quantify 
this finding, we have plotted in Fig.[12]the cumulative distribu- 
tions of two subsamples of stars in the inner (r < 3'.9) and outer 
(r > 3'.9) region, respectively. This limit was chosen to have 
the same number of stars in each subsample. A Kolmogorov- 
Smirnov test indicates that the null hypothesis that the distribu- 
tions are drawn from the same parent population can be rejected 
at > 99.9% confidence level. Similarly, we have plotted the cu- 
mulative distributions of two subsamples of stars with age < 4.6 
Gyr and > 4.6 Gyr. Also in this case, the probability of the null 
hypothesis can be rejected at a level > 99.9%. 

This suggestion of an age gradient among RGB stars is 
strengthened by the detection of a radial gradient in the frac- 
tion of upper- AGB stars in Leo I, which points to a concentra- 
tion of intermedi ate-age populations towards the galaxy centre 
dHeld et al.ll2009l 
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Fig. 12. Ages of the Leo I stars in our sample as a function of the 
elliptical radius (central panel). The upper panel shows the cu- 
mulative radial distributions of stars with ages smaller and larger 
(dotted line) than 4.6 Gyr. The cumulative age distributions of 
stars in the inner (r < 3f 9) and outer (r > 3'9, dotted line) region 
are shown in the right panel (the different age and radial inter- 
vals are colour-coded in the electronic version of the journal). 
Younger stars appear to be more concentrated towards the centre 
of Leo I. 



7. Summary and conclusions 

We have presented spectroscopic measurements of RGB stars in 
the Leo I dSph from observations carried out with the FORS2 
spectrograph at the ESO VLT. We derived radial velocities for 
57 stars with good S/N ratio, 54 of which have been found to 
be Leo I members. Among these, 14 stars are in common with 
previous spectroscopic studies. 

We measured the metallicities of RGB stars in Leo I from the 
equivalent widths of Can triplet lines, usin g the [Fe/H] metallic - 
ity scales of lCarretta & Grattonl d 19971) and lZinn & Westl dl984l) . 
In addition, we derived a new calibration tied to the [M/H] rank- 
ing of Galactic globular clusters, which accounts for the abun- 
dance of both Fe-group and a elements. 

The metallicity distribution (MDF) of Leo I stars is sym- 
metric and very narrow. If we adopt a quadratic calibration 
of Can line strengths against [Fe/H], the mean metallicity is 
[Fe/H] = -1.41 with a m easured dispersion 0.21 dex on the 
ICarretta & Grattonl d 19971) scale, in agreement with previous 
spectroscopic studies. The new [M/H] calibration yields a mean 
value [M/H] = -1.22 with a dispersion 0.20 dex. By subtract- 
ing the measurement errors, we estimated a very low intrin- 
sic metallicity dispersion, <T[m/h] = 0.08, which represents a 
constraint for modelling the chemical evolution of this isolated 
dwarf galaxy. As pointed out by previous studies, this narrow 
MDF is inconsistent with a simple "closed-box" chemical evo- 
lution model, even adopting a very low effective yield to ac- 
count for galactic outflows expelling the metals produced by 
SNe winds. A prompt initial chemical enrichment may explain 
the very small number of extremely metal poor stars (we find 



only 4 stars with [Fe/H] < -2). Together, the two effects can 
explain the small abundance dispersion of Leo I stars, which 
gives the narrowest observed MDF among Local Group dwarf 
galaxies. However, detailed c hemical evolution models (e.g., 
lLanfranchi & Matteuccil [20071) are needed to gain a complete 
picture of the evolution of Leo I. 

Our data for RGB stars also provide an indication of a weak 
radial metallicity gradient in Leo I, of -0.27 dex Kpc -1 . In fact, 
all of our stars with [M/H]> -1.3 are found in the inner region 
(r < 5'). However, by combining our observations with previous 
spectroscopic datasets in the literature, the radial variation be- 
comes insignificant. More observations in the outskirt of Leo I 
with a quality comparable to those presented here, are required 
to definitively establish the presence of an abundance gradient. 

The metallicities of the RGB stars in our Leo I sample have 
been combined with existing photometric data to yield age es- 
timates and an age-metallicity relation. Our age determina tions 
are cons i stent w ith the SFH derived bv lGallart et alj d 19991) and 
iDolphinl (|2002j) from HST photometry. The age-metallicity rela- 
tion of Leo I red giants is quite flat, again suggesting a rapid ini- 
tial enrichment. An increase in metal abundance by ~ 0.2 - 0.3 
dex in the last 5 Gyr is possibly related to the main star-formation 
episode at intermediate ages. Since Leo I only hosts a minor old 
(> 10 Gyr) stellar component, the chemical history of the galaxy 
is not well constrained at early epochs. Its most metal-poor stars 
must have formed out of a medium pre-enriched by a lost gen- 
eration of stars, either before or after the galaxy had started as- 
sembling. 

We have provided the first evidence of a radial variation 
in the ages of red giants in Leo I. Despite the uncertainties in 
age determination, our direct measurement of a radial varia- 
tion of stellar ages seems quite convincing, with a Kormogorov- 
Smirnov test confirming, at a high level of statistical signifi- 
cance, that stars in the inner part of Leo I are on average younger 
than those in the outer regions. This result agrees with the con- 
clusions of a parallel study of inter mediate-age AGB stars in 
Leo I from near-infrared photometry dHeld et al.l l2009t) . In the 
emerging scenario, the first generation of Leo I stars uniformly 
formed throughout this isolated dwarf spheroidal galaxy. The 
bulk of intermediate-age stars originated from an interstellar 
medium, poorly enriched by previous stellar generations mainly 
because of the effects of stellar winds. Younger stellar popula- 
tions preferentially formed in the central regions, from gas some- 
what enriched as seen from the age-metallicity relation in the last 
few Gyr. In this framework, our results on the radial distribution 
of Leo I stellar populations are not in contrast with previous re- 
sults which found no gradien ts. The lack of detection of an age 
gradient bv lHeld et alld2000l) can be explain ed considering tha t 
the mean age of the red-clump stars used by iHeld et all d2000t) 
as tracers of intermediate-age populations, is ~ 5 Gyr, which i s 
older than that of upper AGB stars used by IHeld et all d2009l) . 
As shown in Fig. [12j there are no clear radial variations in the 
age distribution of Leo I stars with ages greater than ~ 5 Gyr. 
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v 
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Av 
km s ~^ 


S/N 


other 


1 
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Notes, v is the heliocentric radial velocity and Av the absolute se mi-difference in radia l velocity of the individual spectra. The last column gives 
the identification of stars in common with lBosler et al .1 d2007l) and lKoch et alj d2007bl) . 
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Table 3. Measurements of metallicity and age for stars in Leo I. 



ID 


ZW 
(A) 


(A) 


W 
(A) 


[Fe/H] zw 


[Fe/H] cu 


[M/H] 


age 
(Gyr) 
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Aage + 
(Gyr) 


1 
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3.71 
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-1.41 
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0.04 


3.86 


-1.48 


-1.34 


-1.15 


5.5 


1.2 
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0.19 


3.19 


-1.77 


-1.63 


-1.43 


9.8 


3.0 


0.9 


9 


5.20 


0.19 


3.85 


-1.48 


-1.35 


-1.15 


5.4 


2.0 


1.3 


10 


2.45 


0.06 


0.82 


-2.17 


-2.14 


-2.04 


6.3 


1.6 


1.9 


11 


4.45 


0.17 


2.85 


-1.89 


-1.75 


-1.55 


6.4 


1.7 


2.0 


12 


4.92 


0.04 


3.57 


-1.62 


-1.48 


-1.28 


6.1 


2.1 


2.7 


13 


5.43 


0.08 


3.81 


-1.50 


-1.37 


-1.17 


2.1 


0.5 


0.8 


14 


5.81 


0.38 


3.98 


-1.42 


-1.28 


-1.10 


2.7 


0.6 


0.8 


15 


4.20 


0.08 


2.84 


-1.89 


-1.76 


-1.55 


9.4 


2.2 


2.4 


16 


5.30 


0.05 


3.76 


-1.53 


-1.39 


-1.19 


5.2 


1.1 


2.1 


17 


4.04 


0.03 


2.60 


-1.96 


-1.83 


-1.63 


4.2 


1.1 


1.6 


18 


6.00 


0.23 


4.16 


-1.32 


-1.19 


-1.01 


2.0 


0.4 


0.3 


19 


5.25 


0.02 


3.72 


-1.55 


-1.41 


-1.21 


5.6 


1.5 


2.6 


20 


4.95 


0.20 


3.26 


-1.74 


-1.60 


-1.40 


1.9 


0.3 


0.4 


21 


5.47 


0.51 


3.51 


-1.64 


-1.50 


-1.30 


3.4 


0.7 


0.9 


22 


5.02 


0.14 


3.59 


-1.61 


-1.47 


-1.27 


11.2 


1.7 


1.4 


23 


6.20 


0.03 


4.34 


-1.22 


-1.10 


-0.92 


4.3 


1.1 


1.9 


24 


5.70 


0.13 


4.17 


-1.32 


-1.19 


-1.00 


3.3 


1.0 


1.5 


25 


5.11 


0.39 


3.67 


-1.57 


-1.43 


-1.23 


4.1 


1.9 


2.0 


26 


4.83 


0.11 


3.50 


-1.64 


-1.51 


-1.30 


5.8 


2.0 


2.1 


27 


5.57 


0.02 


3.83 


-1.49 


-1.36 


-1.16 


5.7 


1.7 


1.8 


28 


5.40 


0.29 


3.82 


-1.50 


-1.36 


-1.17 


4.6 


1.3 


1.5 


29 


4.93 


0.05 


3.12 


-1.80 


-1.66 


-1.45 








30 


2.24 


0.00 


0.75 


-2.17 


-2.14 


-2.05 


10.2 


2.0 


1.0 


31 


5.36 


0.46 


3.72 


-1.55 


-1.41 


-1.21 


3.2 


0.7 


1.4 


32 


4.53 


0.41 


3.05 


-1.82 


-1.68 


-1.48 


6.3 


2.0 


3.0 


33 


4.90 


0.01 


3.10 


-1.80 


-1.66 


-1.46 


8.5 


2.0 


2.2 


34 


5.13 


0.29 


3.29 


-1.73 


-1.59 


-1.39 


1.7 


0.2 


0.3 


35 


5.73 


0.03 


3.91 


-1.45 


-1.32 


-1.13 


4.9 


1.2 


1.7 


36 


5.46 


0.07 


3.57 


-1.61 


-1.48 


-1.27 


2.6 


0.5 


0.8 


37 


5.50 


0.01 


3.71 


-1.55 


-1.41 


-1.22 


10.2 


2.5 


1.9 


38 


5.40 


0.33 


3.68 


-1.57 


-1.43 


-1.23 


7.1 


1.9 


1.8 


39 


5.32 


0.48 


3.82 


-1.50 


-1.36 


-1.17 


2.5 


0.7 


1.3 


40 


5.15 


0.11 


3.66 


-1.57 


-1.44 


-1.24 


7.2 


2.0 


2.7 


41 


5.08 


0.05 


3.21 


-1.76 


-1.62 


-1.42 


5.5 


1.3 


2.0 


42 


6.45 


0.03 


4.34 


-1.22 


-1.10 


-0.92 








43 


6.12 


0.12 


4.45 


-1.16 


-1.04 


-0.87 


2.8 


0.8 


1.2 


44 


5.28 


0.21 


3.44 


-1.67 


-1.53 


-1.33 


2.6 


0.7 


0.7 


45 


2.94 


0.23 


1.52 


-2.16 


-2.07 


-1.92 


4.9 


1.5 


1.9 


46 


6.13 


0.07 


4.33 


-1.23 


-1.10 


-0.93 


3.8 


1.1 


1.3 


47 


5.86 


0.44 


4.07 


-1.37 


-1.24 


-1.05 


1.5 


0.1 


0.2 


48 


6.21 


0.21 


4.51 


-1.12 


-1.00 


-0.84 


1.7 


0.2 


0.2 


49 


3.26 


0.26 


1.41 


-2.16 


-2.09 


-1.95 


8.0 


1.9 


1.9 


50 


5.81 


0.04 


3.65 


-1.58 


-1.44 


-1.24 








51 


4.96 


0.06 


3.09 


-1.81 


-1.67 


-1.47 


7.8 


1.3 


1.0 


52 


5.79 


0.38 


3.89 


-1.46 


-1.33 


-1.13 


4.0 


1.0 


1.1 


53 


5.67 


0.11 


3.96 


-1.43 


-1.29 


-1.10 


4.6 


1.4 


2.4 


54 


6.17 


0.22 


4.55 


-1.09 


-0.98 


-0.81 


1.9 


0.3 


0.7 



Notes. e-£.w is the absolute semi-difference of the equivalent widths measured on the individual spectra. The listed metallicity values were 
calculated from W using the quadratic calibration. The last two columns give the lower and upper confidence intervals of our age estimates (see 
text for details). 



